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The formation of intermolecular bonds in Cgo has been investigated in detail at pressures below 2.2 GPa and 
up to 750 K. FuUerene samples were heated in a temperature gradient to obtain data on the formation of 
dimers and low-dimensional polymers along isobars. Intermolecular bonding was analyzed ex situ by 
Raman scattering, using both intramolecular modes and intermolecular stretching modes. 
Semi-quantitative reaction maps are given for the formation of dimers and chains. The activation energy for 
dimer formation decreases by 0.2 meV pm~' when intermolecular distances decrease and dimer formation 
is noticeably affected by the rotational state of molecules. Above 400-450 K larger oligomers are formed; 
below 1.4 GPa most of these are disordered, with small domains of linear chains, but above this the 
appearance of stretching modes indicates the existence of ordered one-dimensional polymers. At the highest 
pressures and temperatures two-dimensional polymers are also observed. 

The physical properties of Buckminsterfullerene, Cgo. can be tuned over wide ranges by high pressure 
treatment. Although the Ceo molecule is much less compressible than diamond'-^, the weak intermolecular 
interactions make solid bulk Ceo three times more compressible ("softer") than graphite'*. However, 
irradiation of solid Cgo by light"" or particles, or compression of the material'', changes its properties dramatically. 
In both cases, covalent intermolecular bonds may form by a [2 + 2] cycloaddition reaction in which reactive C = C 
sp*" ("graphite-type") intramolecular bonds on nearest neighbours break up to form two parallel intermolecular 
sp"" ("diamond-type") C-C bonds. The polymerization reaction changes the mechanical and electrical properties 
and high pressure treated C^o has been reported to have band gaps from^ almost zero to" 2.35 eV, while its bulk 
modulus and hardness approach, or possibly even exceed', those of diamond. 

The possibility to tune the properties of Ceo over wide ranges in this way has attracted much interest. Pressure- 
induced polymerization has been extensively studied*"" by treating Ceo at different pressure-temperature con- 
ditions and measuring its structural and vibrational properties both in situ and ex situ. Some investigations have 
involved over a hundred individual high pressure experiments" '^ and the results have been used to draw reaction 
maps for the products obtained under different conditions" The experiments show that compression at room 
temperature produces random intermolecular bonds, initially at a low rate. However, near 10 GPa a more rapid 
reaction is reported^" followed by a sudden lattice coUapse^"'^' near 20 GPa, probably due to massive three- 
dimensional (3D) intermolecular bond formation. The original lattice can be stabilized to higher pressures by 
inserting metal ions^^ or solvents^' as spacers to prevent the reaction. At elevated temperatures and modest 
pressures (1-5 GPa), pure Ceo forms several well defined covalently bonded structures'"'"''" ''^ ("polymers"). 
The first step is the formation of (C6o)2 dimers, but with increasing pressure and/or temperature larger oligomers 
grow. Below 10 GPa and 1000 K, three long-range ordered structures are found"''" "''": an orthorhombic phase 
containing parallel linear molecular chains (one-dimensional polymer, ID), a pseudo-tetragonal phase, where 
such chains are linked by bonds in the perpendicular direction, and a rhombohedral phase where all molecules in 
the original (111) planes of the face-centered cubic (fee) Ceo lattice are covalently bonded. The last two both 
contain two-dimensional (2D) polymer planes. Once formed, polymers are metastable at room temperature but 
break down to monomers on heating at atmospheric pressure; they are stable to higher pressures^* than mono- 
meric Ceo but may transform into 3D polymers under extreme conditions'''^" '^'^'"^^. 

Although the general features of this polymerization process are known, there is surprisingly little agreement 
on details. The main reasons are the random disorder inherent in the process, which complicates the structural 
analysis, and uncertainties in the pressures and temperatures. Pressure is often calculated from applied load and 
temperature is measured using thermocouples. Both methods give a reasonable accuracy, but in small pressure 
cells small changes in the dimensions or materials properties may give significant changes in pressure, heating 
changes the pressure by softening the pressure transmitting medium or by thermal expansion, and large tem- 
perature gradients up to IC to 10^ K m"' increase uncertainties in temperature measurements. Polymerization is 
also sensitive to anisotropic strain which promotes the reaction by locally bringing molecules into contact. 
Furthermore, X-ray diffraction or Raman scattering results are often presented in a non-quantitative way. 
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Figure 1 | Raman scattering data, (a) Raw Raman data for four samples after pressure treatment. As indicated by the curve labels, spectra are shown 
both for the pristine material, a sample dominated by dimers, a sample dominated by chains, and a sample containing both chains and a high fraction of 
2D tetragonal Cgo ("planes"). Top two curves have been magnified and shifted vertically for improved clarity but the original background has been 
preserve to indicate data quality, (b) Fitted components of the Ag(2) mode for the same four samples as in (a), after removal of baseline. The peaks have 
been reseated to approximately the same total amplitudes. The fitted curves are labeled pristine, dimer, chain, fork and plane, indicating signals 
originating from free molecules and from molecules covalently bonded to 1, 2, 3 and 4 neighbours, respectively. 



identifying only the majority phase. Creating a reliable reaction dia- 
gram from the results of experiments using samples of different 
purity, different high pressure devices, reaction pressures, pressure 
media, temperatures and reaction times, analyzed by different meth- 
ods and treated at widely different pressure-temperature coordi- 
nates, is thus a challenging task. 

Recently, a simple method was developed^" to map reaction dia- 
grams by heating a long sample in a temperature gradient under 
pressure, then quenching to ambient conditions for analysis. Each 
experiment investigates an isobar in the pressure-temperature plane 
and a small number of experiments provide a complete two- 
dimensional view. For metastable phases each experiment provides 
a semi-permanent record which can be studied ex situ, and interest- 
ing regions can be re-analyzed in detail if desired. Although the 
absolute temperature and pressure accuracy is limited by the same 
factors as in other experiments, the pressure is constant along each 
sample and relative temperatures are extremely well defined. Here, 
this method has been applied to the reaction diagram of Ceo- The 
results have been analyzed using Raman scattering, focusing on the 
intermolecular bonding rather than on lattice structure, and a semi- 
quantitative map of the reaction products at any point in the 
pressure-temperature plane has been obtained. The data verify that 
the rotational, and probably also the orientational, states of the mole- 
cules influence the intermolecular bonding reaction. 

Experimental results 

This report is based on data from 16 samples, each treated for one 
hour at pressures from 0.7 to 2.2 GPa. Except at the lowest tempera- 
tures this should be long enough for the reaction to reach equilib- 
rium. The heater power was usually chosen to give a temperature 
range from about 375 to 550 K, but three samples were run with 
maximum temperatures near 750 K to study the onset of two- 
dimensional polymerization. As described in the Methods section, 
samples were well characterized and great care was taken to avoid 
photo-induced effects before and after high pressure treatment. 

Raman spectra. The presence of intermolecular bonds has been 
identified using NIR Raman scattering. Raman spectra were 
measured at 20-150 points on each sample and the corresponding 



treatment temperatures calculated from the measured mechanical 
position. Including tests on pristine samples, 841 Raman spectra 
were measured, most of these using 830 nm laser excitation. 

The number of intermolecular bonds to each Ceo molecule is 
determined from the shift of the Ag(2) intramolecular Raman 
mode'^", but this method provides no information about lattice 
structure. However, it was also possible to observe the longitudinal 
lattice mode of linear Ceo chains near 117 cm"' and the dimer 
stretching mode at 96 cm"'. Fig. la shows Raman data for four 
samples: the pristine material before pressurization, a sample con- 
taining mainly dimers, another containing a high concentration of 
the "linear chain" phase and a fourth containing both chains and 2D 
polymer. For the pristine material all ten Raman modes allowed by 
symmetry, two Ag modes (at 495 and 1469 cm"') and eight Hg 
modes, are easily identified, and because isotopes and defects break 
the local symmetry some further peaks may be observed in the noise 
(not visible on this scale). Although ordered dimer lattices have been 
predicted^''" and may have been observed"", most experiments 
show""'™''^ that dimers form at random in the lattice. The resulting 
disorder breaks the lattice symmetry creating a very rich Raman 
spectrum. When ordered polymer lattices form, the symmetry again 
becomes higher and the number of peaks decreases. The evolution of 
the Raman spectrum with polymerization has been well investi- 
gated'^'"*'''^'''^""'. Features of particular importance for the present 
work are: i) the Ag(2) mode at 1469 cm"' shifts in a predictable 
way with the number of sp"" covalent intermolecular bonds: for each 
cycloaddition link formed, the peak shifts by approximately 
— 5 cm"', ii) The Ag(l) radial breathing mode at 495 cm"' shifts 
by about —6 cm"' on dimerization. iii) Polymerization also gives 
rise to new modes; in particular, both dimers and chains show low 
frequency longitudinal stretching modes, and iv) the formation of 
intermolecular bonds causes new modes in the range 900- 
1000 cm"', the so called "polymer modes" associated with'' vibra- 
tions in the four-membered ring of the intermolecular bonds. The 
first three features are used here to identify the polymeric structures 
after treatment, and polymer modes were also observed: lines near 
950 and 975 cm"' appear when dimers form, and with increasing 
polymerization these lines broaden and merge into a band with 
several peaks, extending from about 940 to 1010 cm"'. With NIR 
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excitation these lines were weak and no systematic study was possible 
but with visible light excitation their intensity is similar to that of the 
Ag(l) or low Hg modes and clear structures are seen. 

The Ag(2) modes from the spectra in Fig. la are analyzed in detail 
in Fig. lb. Voigt functions with peak positions near 1469, 1464, 1459, 
1454 and 1448 cm"', implying covalent intermolecular bonds to 0-4 
nearest neighbours, respectively, have been fitted to the data. For 
simplicity these peaks will be referred to as originating from mono- 
mers, dimers, chains, forked (branched) chains and planes, respect- 
ively. All measured spectra were analyzed in this way and in most 
cases the signal-to-noise ratio was sufficient to give stable values for 
the relative peak areas. The peak near 1448 cm"', associated with the 
tetragonal 2D polymer, was only found in samples treated at the 
highest pressures and temperatures and was usually associated with 
a line near 1432 cm"', believed'"' to be a sub-peak of the Hg(7) mode. 
A small peak near 1454 cm"' often significantly improved the fit but 
does not correspond to any known lattice structure. Its position 
suggests bonding to three molecular neighbours, i.e. a structure with 
branched chains or oligomers. The relative fitted areas of the Ag(2) 
sub-peaks are used as semi-quantitative measures of the numerical 
fractions of molecules with different numbers of intermolecular 
bonds. Accurate Raman cross sections for different Ceo structures 
and excitation wavelengths are not known, since pure structural 
phases are rare, and it is thus not possible to deduce the true fractions 
of the structures observed. However, where measurements were car- 
ried out using all three excitations lasers the results usually agreed to 
within the experimental scatter. Figure la shows that the sensitivity 
of the detector drops off rapidly with increasing wavenumber near 
the Ag(2) mode, which might overestimate the intensity from poly- 
merized Ceo relative to monomers. Compared to other possible 
errors this is no large effect and no correction for this has been 
applied. 

For each experiment, the fitted relative intensities from different 
molecular states were plotted as functions of treatment temperature 
as shown in Fig. 2. At low temperatures all samples contained mainly 
Ceo monomers, but a weak dimer line was usually observed. 
Although relative intensities at the single percent level are difficult 
to quantify by fitting, the spectra always showed indications of other 
dimer peaks in the noise when the 1464 cm"' peak reached this level. 
In the intermediate pressure range, 1.1-1.8 GPa, the relative intens- 
ities evolved in an interesting multi-step sequence. At low temper- 
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Figure 2 | Relative intensities of Ag(2) sub-modes versus treatment 
temperature. The relative intensities of the five Ag(2) sub-modes for a 
sample treated at 1.4 GPa, plotted as functions of treatment temperature. 
Signals originate from monomers and from Cgo molecules covalently 
bonded to one ("dimer"), two ("chain"), three ("fork") and four 
("plane") nearest neighbours. 



ature the intensity from dimers was practically constant. At some 
threshold temperature this intensity started to grow and a small 
constant intensity from chains appeared. When the area of the dimer 
peak reached 20-35 percent of the total Ag(2) area, a sudden step 
increase in dimer intensity by about a factor of two occurred. This 
step was very sharp, with a width of often less than 1 K, from 1.2 to 
1.6 GPa. Above the step the dimer intensity was constant or 
increased slowly with temperature, the intensity from monomers 
decreased rapidly and the intensity from chains increased approxi- 
mately linearly with increasing temperature. Usually the dimer 
intensity showed a plateau or a rounded maximum, 20-40 K wide, 
before dropping off at the same rate as the chain intensity increased. 
At the highest temperatures the NIR Raman measurements showed 
constant, small relative intensities from monomers and dimers, a 
dominating contribution from chains, and for the higher pressures 
significant contributions from "branched chains" and 2D tetragonal 
polymers. No signals were found from the rhombohedral phase. 
Considering that polymerization is a random process the scatter in 
the data was surprisingly small. However, in most samples grains 
were found which differed strongly from their neighbours, possibly 
because of local strain during treatment. 

Reaction map for dimers. The relative intensities of the Raman 
peaks may be plotted as functions of temperature and pressure to 
create semi-quantitative reaction maps for each molecular state. 
Although the scatter in the data (Fig. 2) suggests that the map will 
be rather crude, it may still be significantly more detailed than earlier 
maps""'^ because points can be spaced arbitrarily close in 
temperature on each isobar. Figure 3a shows a map of the relative 
intensity of the 1464 cm"' dimer line. The diagram was drawn by 
connecting the points in graphs of intensity versus temperature 
(Fig. 2) by straight lines and noting where these lines cross selected 
percentage levels. To minimize artifacts the curves were not 
statistically smoothed and the experimental scatter is intentionally 
preserved. Most contours in Fig. 3a were interpolated between 
experimental points using cubic spline functions, manually 
modified when obvious artifacts, such as crossing curves, occurred. 
The highest intensities are found only in small areas at a few 
pressures, and for these the contours have been drawn manually. 
The scatter in the data also creates small hills or valleys, again 
drawn manually. Data were extrapolated by 50-100 MPa for the 
final map to cover the range 0.65 to 2.3 GPa. 

Because dimers are often considered as precursors to larger oligo- 
mers most earlier reaction maps"""" do not show under what con- 
ditions they form, and even when dimerization is discussed such 
information is qualitative only"'^^"" '". In contrast. Fig. 3a clearly 
shows that dimers are the major reaction product when Ceo is 
annealed for one hour under certain pressure-temperature condi- 
tions. Both formation and breakdown of intermolecular bonds are 
believed to be thermally activated processes and the observed reac- 
tion diagram is in good general agreement with this assumption. The 
dimer Raman intensity increases with reaction temperature, until the 
dimers start coalescing into larger oligomers and the intensity drops. 
However, a surprisingly large intensity (10-30 percent) remains at 
high pressure and temperature, where larger oligomers should 
dominate. 

The Ag(2) line originates from vibrations along the C = C bonds on 
individual molecules and a shift of this line might be caused by other 
mechanisms than dimerization. Another way to detect dimers is 
through their Raman active stretching mode, low frequency longit- 
udinal vibrations of the heavy molecules along the intermolecular 
bond. The dimer spectrum in Fig. la was selected to show a strong 
stretching mode at low wavenumber, 96 cm"'. The intensity was 
quantified by calculating the ratio between the intensities of this 
mode and of the nearby radial breathing Ag(l) mode. Because the 
latter shifts from 495 cm"' for the monomer to 489 cm"' for dimers 
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Figure 3 | Reaction map for dimer formation, (a) Reaction map obtained from the fitted intensities of the molecular Ag(2) mode. Colour shows the fitted 
relative intensity of the 1464 cm"' component of the Ag(2) mode as a function of pressure p and temperature T. (b) Reaction map obtained by plotting the 
ratio Idimer/lAg(i) between the intensities of the longitudinal stretching mode at 96 cm"' and the total intensity of the Ag(l) breathing mode at 489- 
495 cm"'. The dashed lines indicate the transition between the low temperature phase with orientationally ordered molecules and the high temperature 
phase with quasi-free molecular rotation, extrapolated from low temperatures using literature data*'". 



and linear chains, the sum of the intensities of these two sub-modes 
was used. A weak third line near 485 cm"', tentatively identified"" as 
Gg(l), usually had to be included in the fit but was not included in 
calculating the ratio. Figure 3b shows a map of the intensity ratio 
Idimer/lAgi as a function of temperature and pressure, plotted using 
the same procedure as for the Ag(2) line; note that a different numer- 
ical scale has been used because of the low intensities. The two 
reaction maps in Fig. 3 are virtually identical, except that Fig. 3b 
contains areas at both high and low temperature where no signals 
were detected. 

Figure 3 shows that dimers are produced in high yield in a region 
extending from 430 K near 1 GPa to 400 K near 2 GPa. At lower 
temperatures few dimers form, and at higher temperatures larger 
oligomers grow. The results from inter- and intramolecular modes 
agree that the highest intensities, above 80%, are obtained at rather 
low pressures, 1.15 ± 0.05 GPa, probably because the reaction rate, 
and thus the formation of larger oligomers, increases with pressure. 
The peak yield is comparable to that found in matrix-assisted dimer 
synthesis'"'. Five features merit special comments: 



i) The dashed line indicates the phase line for the rotational 
transition found near*' 260 K at atmospheric pressure, extra- 
polated using literature data'^'*^. Surprisingly high intensities 
were found near this line for both molecular and stretching 
modes even at quite low temperatures, and the data suggest 
the existence of a low pressure "tail" of the dimerization reac- 
tion close to this transformation. (Note that a small shift 
( + 100 MPa) of the position of the line would bring it into 
excellent agreement with areas of local maximum intensity 
from 0.8 to 1.2 GPa; the normal phase line curvature should 
modify the linear extrapolation in this direction.) The resolu- 
tion of earlier experiments has been too low to detect effects of 
molecular rotation and orientation on the rate of intermol- 
ecular bonding. However, the dimerization enhancement is 
unexpected. Although the reaction rate increases with tem- 
perature the rotational transformation is associated with a 
lattice expansion which should decrease the reaction rate. 
The high dimerization rate at low temperatures below 
1 GPa rather suggests that transformation-induced dynamic 
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disorder locally brings molecules into close contact, enhan- 
cing the reaction rate. 



large, Figs 3a and 3b agree that the temperatures for both 
dimerization onset and maximum dimer intensity change 
with pressure at a rate of ( — 28 ± 5) K GPa"'. The more 
interesting relation between dimerization temperature and 
intermolecular distance is presented in Fig. 4. Defining dimer- 
ization onset temperatures from the sharp increase in dimer 
intensity discussed above and using intermolecular distances 
in Cgo from an internally consistent set of compression 
data"'*^'''^ and a linear coefficient of thermal expansion'"' of 
1.52* 10"^ K"' gives a well-defined straight line valid for both 
rotating and non-rotating molecules. Here the volume change 
at the rotational transition has been taken into account, but it 
should be remembered that the definition of "dimerization 
temperature" is arbitrary and that the lattice spacing is over- 
estimated, since the sample at this point contains 20-35 per- 
cent dimers. Extrapolation shows that at 293 K the onset of 
rapid dimerization should occur at a volume compression of 
20 percent, in excellent agreement with the report of signifi- 
cant diffraction line broadening by Haines and Leger™ and an 
anomaly in the data of Horikawa et al.*^ at this compression. 
The slope 2.45 K pm"' indicates that when fullerene 
molecules are pushed together the bond formation energy 
threshold decreases by about 0.2 meV pm"' in the range 
investigated. As expected from the nearly isothermal dimer- 
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Figure 4 | Dimerization onset versus intermolecular distance. The 

characteristic temperature for the sudden onset of intense dimerization, as 
defined in the text, plotted as a function of estimated nearest-neighbour 
spacing under the corresponding pressure-temperature conditions. The 
three points at the largest spacing corresponds to a face-centered lattice 
with almost free molecular rotation, all other to the low-temperature 
simple cubic lattice with partial orientational order. 



ization behavior this is small compared to the activation 
energy, for which calculations''^''"' show values of 3-4 eV while 
experimental data range from'" 0.4 to'"' 1.4 eV. 
A second anomaly may exist near 2 GPa, where the dimer 
intensity drops off rapidly. This may be correlated with the 
transition'^ into an orientationally ordered lattice, which from 
a linear extrapolation of low- temperature data' "" should occur 
near 1.5 GPa at 350 K and 2 GPa at 500 K. Again, the normal 
phase line curvature should push these values to higher pres- 
sures. The absence of molecular ratcheting in the ordered 
lattice might delay initial dimerization to higher temperature. 
Once intermolecular bonding starts, however, molecular ori- 
entation might facilitate growth of larger oligomers, leading to 
a decrease in maximum dimer intensity, in agreement with 
the data. 

Finally, above 450 K the Ag(2) mode (Fig. 3a) indicates sig- 
nificant dimer intensities but there is practically no signal 
from stretching modes (Fig. 3b). There is no obvious explana- 
tion for this difference, but a strong 1464 cm"' Ag(2) line is 
often observed also in tetragonaP* (2D) Ceo- The reason is 
unknown; possibly, surface reconstruction adds a high density 
of terminal molecules attached to the edges of planes by 
cycloaddition bonds. 

Reaction map for fullerene chains. Using the same procedure, 
similar maps were drawn for Ceo molecules covalently bonded to 
two neighbours, usually associated with an orthorhombic lattice 
containing linear chains'" ''^. The map in Fig. 5a shows the relative 
intensity of the 1459 cm"' Ag(2) sub-line, extrapolated to 0.65 and 
2.3 GPa as before. Three experiments were extended to about 750 K. 
The selected pressures were 1.15 GPa, to check the position of the 
reported equilibrium phase line^° between Cgo monomer and 
polymer, 1.3 GPa, because early runs indicated that 2D structures 
begins to form here, and 1.9 GPa, to study the evolution of 2D 
polymers. White areas denote either that no signals were detected 
(low temperatures) or that no measurements were made (high 
temperature) and the dashed line is again the (extrapolated) 
rotational transition boundary of molecular Ceo- A comparison 
with Fig. 3a shows that the decreasing dimer Raman intensity is 
matched by an increase in intensities from chains, rising steeply at 
temperatures between 400 and 450 K. Near 1 GPa, 60-85 percent of 
the Ag(2) intensity originates from chain molecules above 450 K and 
we showed earlier^' that under these conditions well ordered ID 
polymers formed in large Ceo single crystals. Above 1.6 GPa the 
chain intensity is lower and with increasing temperature it rapidly 
drops to below 50 percent, suggesting that chains link up to form a 
2D polymer. These results are in excellent agreement with literature 
jj^(jjii-i7,3o.3335.38,5i^ cspeclally the work of Davydov et al." 

Like dimers, chains can also be identified from a longitudinal 
lattice mode (near 117 cm"') and Fig. 5b shows the corresponding 
map, drawn in the same way as for the dimer stretching mode 
(Fig. 3b). It is immediately obvious that Figs 5a and 5b are very 
different. In the high intensity area of Fig. 5a no stretching mode 
intensity is observed. The stretching mode intensity instead peaks 
near 2 GPa and 600 K. (Note that the figure shows the ratio relative 
to the radial breathing Ag( 1 ) mode which rapidly loses intensity with 
an increasing number of intermolecular bonds.) A reasonable 
explanation for the difference between the figures is found by real- 
izing that the lattice stretching mode is only excited if sufficiently 
long, straight chains exist, while the intramolecular mode is the same 
for all molecules bonded to two neighbours. Combining the present 
data with earlier results"'^''^' suggests that below 1.4 GPa dimers 
coalesce to form small (200 A) ordered domains^' of short chains 
in a matrfac of disordered oligomers like non-linear, randomly zig- 
zagging chains. A significant (5-15 percent) intensity from an Ag(2) 
sub-line near 1454 cm"' suggests that such oligomers form branched 



ii) In the low-temperature simple cubic phase the onset of rapid 
dimerization is surprisingly sudden between 1.1 and 1.8 GPa. iv) 
The intensity from dimers jumps from 20-30 to 50-60 per- 
cent for a temperature increase of less than 1 K (Fig. 2). This is 
particularly surprising considering the drift in treatment tem- 
perature with time (see section Methods). Below this jump the 
intensity from the Ag(2) dimer line increases approximately 
linearly with temperature, while above it is almost constant. 
The rapid increase might be caused by positive feedback. 
Dimerization locally reduces the lattice volume and the result- 
ing dynamic disorder may allow other molecules to take up 
positions favourable for further dimer formation. 

iii) Pressure forces molecules together, decreasing the dimeriza- 
tion onset temperature. Although the experimental scatter is v) 
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Figure 5 | Reaction map for chain formation, (a) Reaction map obtained from the fitted intensities of the molecular Ag(2) mode. Colour shows the fitted 
relative intensity of the 1459 cm"' component of the Ag(2) mode as a function of pressure p and temperature T. (b) Reaction map obtained by plotting the 
ratio I(-i,ain/lAg{i) between the intensities of the longitudinal chain stretching mode at 117 cm"' and the total intensity of the Ag(l) breathing mode at 489- 
495 cm"'. The dashed lines in these figures are identical to those shown in Fig. 3 while the full lines indicate the equilibrium between polymeric and 
monomeric Ceo phases as given by Bashkin et al?° The latter lines have been extrapolated linearly to low temperatures. 



structures. Both the dimer and chain "phases" are thus strongly 
disordered below 1.4 GPa and large crystallites of ordered ortho- 
rhombic (chain) phase should only be expected at higher pressures 
where the stretching mode is observed (Fig. 5b). 

The full lines in Fig. 5 are extrapolations from high temperature of 
the equilibrium line between monomers and polymers reported by 
Bashkin et alJ^ Figure 5a shows excellent agreement with this line 
near 750 K; near the line the intensity drops rapidly. (The parallel 
arrangement of the contour lines is a reasonable guess.) Near 500 K, 
Bashkin et al. found about 50% polymerization at 0.8 GPa; Fig. 5a 
shows a higher polymer fraction, but considering hysteresis effects 
the agreement is quite good. Again, effects of the rotational state 
might be present. Below the rotational phase line (dashed) the poly- 
merization onset shifts with pressure at about the same rate as for 
dimerization, but above the line the slope is an order of magnitude 
larger, 3-400 K GPa"', indicating an intermolecular distance too 
large for efficient bond formation. 

Two-dimensional polymers. Above 1.2 GPa a further Ag(2) sub- 
mode can be detected near 1448 cm"', indicating that each 
molecule has formed bonds to four neighbours. This line is usually 
associated with the pseudo-tetragonal"'"'^'''^^"^"' 2D phase. The NIR 



signal-to-noise ratio rapidly degraded with increasing treatment 
temperature and pressure, but control experiments using visible 
lasers gave very similar results. Fig. 6 shows the relative intensity of 
the 1448 cm"' line as a function of treatment temperature and 
pressure. At and below 1.2 GPa this mode was just detectable at a 
few points, but at 1.3 GPa weak lines were clearly observed above 
520 K. At and above 1.4 GPa the intensity increased rapidly and 
linearly with temperature, with the onset temperature gradually 
falling from 480 K to near 400 K. The observation of such low 
onset temperatures and pressures for the formation of the 2D 
phase was unexpected; most earlier studies observe this phase only 
under much more severe conditions. One reason might be the long 
reaction time used; the tetragonal phase may grow from the linear 
chain phase^^ and one hour should be long enough for the material to 
evolve from monomers into dimers, linear chains and finally 2D 
planes. Although Fig. 6 indicates maximum relative intensities 
below 50 percent, the true 2D polymer fraction may be higher 
because the Raman response was found to weaken with increasing 
polymerization. The results agree well with our earlier studies: X-ray 
diffraction indicated^^ that a single crystal treated for four hours 
under hydrostatic conditions at 700 K and 2 GPa was almost 
completely converted into 2 D polymers. The tetragonal phase was 
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Figure 6 | Relative intensity of the Ag(2) Raman mode from 2D polymer. 

The relative intensity of the Ag(2) sub-mode associated with molecules 
covalently bonded to four nearest neighbours, believed to indicate the 
presence of the tetragonal 2D polymer of Ceo) as a function of temperature 
at several pressures. Three points deviating very far from the general 
behaviour have been removed for clarity. 

the majority phase, but in contrast to the present samples a small 
amount of rhombohedral phase was also foimd. Raman studies on 
the polymerized single crystals^^ also showed signals from dimers 
and linear chains. In particular, the dimer Ag(2) peak at 1464 cm"^ 
was quite strong, in agreement with the data shown in Fig. 3a above. 

Discussion 

In contrast to most earlier studies of the polymerization of Ceo this 
work focuses on the early stages of intermolecular bond formation 
rather than on the lattice structure of the resulting material. In par- 
ticular, it was shown that (050)2 dimers are a stable, dominating 
structural feature in a certain temperature interval at all pressures 
above 800 MPa, and that the rotational transition is associated with 
an enhanced dimer formation rate, possibly due to dynamic disorder. 
The results strongly emphasize the inherent randomness of the poly- 
merization process, implying that no pure structural phases are 
obtained in any pressure- temperature region studied here. In par- 
ticular, at "low" pressure (below about 1.5 GPa) polymerization of 
Ceo always results in very strong structural disorder. Above 450 K 
randomly oriented dimers grow or merge by rapid formation of 
intermolecular bonds, producing a disordered structure consisting 
of oligomers with a random structure and little long-range order. In 
this range most molecules are connected to two or three neighbours, 
probably forming networks of branching, non-linear chains, and the 
absence of lattice stretching modes indicates that only few and/or 
small domains of ordered orthorhombic polymers exist. These fea- 
tures may be invisible in X-ray studies, which "see" only long-range 
order. Surprisingly, two-dimensional polymers are detected already 
near 1.3 GPa and 500 K, conditions much more gentle than com- 
monly believed necessary for the formation of such structures. The 
detailed information presented here thus gives new possibilities to 
improve our understanding of thermally induced polymerization of 
Ceo; however, much remains to learn, particularly considering the 
fact that irradiation with light is known to further affect both thermal 
polymerization^^ under pressure and thermal depolymerization^^. 

Methods 

Sample material. The material investigated was pre- sublimed, nominally 99.98 
percent pure Cgo obtained from Term USA, CA. Before use the material was slowly 
heated to 200 "C and kept at this temperature for 24-48 h in a dynamic vacuum {10~* 
mbar) to remove atmospheric gases and water, and to break up any photopolymers 



formed during storage. The samples were removed from the vacuum chamber in a 
glove-box and stored under dark conditions in argon untU use. 

High pressure cell. All samples were treated in a piston -cylinder device, 45 mm in 
inside diameter and similar to that described by Andersson et al.^^ In each cell, the 
sample was distributed in the form of powder on a 4 mm wide strip cut from 0.3 mm 
thick copper foil, placed horizontally in the center plane of the Teflon"^ pressure cell. 
The amount of material was chosen to form a 200 |im thick layer after compression. 
The copper substrate was heated at one end using a Kanthal^*^ wire heater on 
pyrophyllite insulation. The useful sample length was usually at least 25 mm. Six 
Constantan wires, 0.25 mm in diameter, were placed on the copper surface with wire 
positions as close as possible to being equally spaced and perpendicular to the 
substrate axis. A 35 |,im Cu foil strip was placed on top of this assembly before a layer 
of talc powder was added to fill out voids and minimize non-hydrostatic pressure. 
Although the intended function of this foU was to minimize mixing of the sample with 
the talc powder, it also proved to be a very convenient aid in determining the final 
positions of the thermocouples, since the Constantan wires usually left clear 
impressions in its surface. The area close to the heater was packed with (heat- 
resistant) talc powder, usually mixed with glass beads to improve thermal insulation^^. 
Finally, the top half of the cell was filled with Teflon, to create a symmetrical cell 
structure, and a top layer of talc powder. To measure the temperature at six points, a 
copper wire was soldered to the cold end of the substrate and used with the six 
Constantan wires to form six thermocouples. Thermocouples made from identical 
materials were calibrated against a calibrated diode sensor (Lake Shore Cryotronics) 
up to 175 C and a standard Pt-100 sensor above. Pressure was estimated from the 
load, using a calibration function derived from measurements in a fluid medium 
using a Manganin wire calibrated against the melting point of mercury. A photograph 
of a typical pressure cell before sample loading is presented as Figure Si in 
Supplementary Material. 

High pressure treatment. All samples were nominally treated for one hour at the 
final pressure- temperature conditions; because the rate of approach to equilibrium 
varied, the effective treatment period was a few minutes longer in some experiments. 
After one hour at the reaction conditions the heater power was switched off, and after 
rapid cooling to room temperature the pressure was slowly decreased to zero. Because 
the final reaction products of Ceo depend more on temperature than pressure, the low 
pressure experiments were carried out by first increasing the pressure to 0.5 GPa, 
then heating to the final temperature, and finally increasing pressure to the final value. 
The reaction then proceeds directly from the monomer phase to the final structure, 
and earlier work on single crystals^^'^^ and thin fUms^" showed that this gave better 
phase purity. However, increasing pressure at high temperature increases the risk for 
broken connections and for all pressures above 1.15 GPa the traditional path, with 
pressurization followed by heating, was used. At 1 GPa two experiments were made, 
one using each method. The results were even quantitatively very similar, except that 
the results for the sample compressed before heating were consistently shifted to 
about 10 K higher temperature compared to the sample heated before compression. 
The heater power was 30-80 W, depending on heater design and the reaction 
conditions chosen. This was enough to give a significant temperature increase in the 
whole vessel with time. The indicated temperature at a single thermocouple could 
easily be kept constant within better than 0.5 K by slowly decreasing the heater power, 
but this led to a slow cooling of the hot part of the sample and a slow heating of the 
cooler part. Neither effect was very large; the maximum change noted over one hour 
was about 5 K. However, the effect complicated the definition of the exact reaction 
temperature. Noting that the polymerization reaction usually runs to near completion 
in a few minutes^^, and that the polymerization is not easUy reversible^^, the solution 
was to define the reaction temperature at each point as the average value during the 15 
minutes when the measured temperature was the highest. Usually this meant the first 
15 minutes for the hot part of the sample and the last 15 minutes for the cool part. The 
actual temperature at each point was then calculated from the position relative to the 
nearest thermocouples, usually interpolated using a second order polynomial. 

Analysis. To avoid post-treatment photo-polymerization all pressure cells were 
opened under red light in a photographic darkroom. Most cells could be cut open and 
disassembled without much sample damage. The Constantan wires conveniently 
divided each pressed sample into five rectangular platelets which could be carefully 
removed individually. Each platelet was then fixed to a glass plate coated with glue 
(with the heated side up) and stored in a dark box. In some cases platelets broke and 
pieces were lost, but for most samples a reasonably complete record of the reaction 
process could be obtained. The reaction conditions at each point on the platelets could 
then be found from the measured temperatures and the position. The estimated 
temperature gradients varied between 3 and 50 K mm~\ and with a simple 
micrometer stage the highest theoretical temperature resolution was thus of the order 
of 10 mK. AU samples were analyzed using a Renishaw 1000 filter-based Raman 
spectrometer with three excitation lasers (wavelengths 514, 632 and 830 nm). All 
measurements presented in the paper were carried out using 830 nm excitation, 
because the energy of this laser is not high enough to cause photo-polymerization. 
The 514 nm laser gave very clear spectra for material polymerized at high 
temperature, but for monomer samples a significant degree of additional 
oligomerization could be observed as soon as the exposure time was long enough to 
provide a reasonable signal-to-noise ratio. The 632 nm laser gave very strong 
fluorescence signals for monomeric Ceo but produced good signals for well 
polymerized material. The 830 nm NIR laser gave excellent spectra for pristine Cgo as 
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well as for most polymerized phases, but for two-dimensional polymers or partially 
decomposed ("amorphized") Cso it gave weak signals with a very strong background. 
A small thermal drift was observed in the laser wavelength, causing spectra to shift by 
about ± 1 cm"' during a day, but since peak positions are well known this was not a 
problem. The sensitivity of the detector also decreased significantly at long 
wavelengths, and with NIR excitation both signal and background fell rapidly for 
Raman shifts larger than 600 cm"'. Spectra were still good enough for accurate 
analysis even at the "pentagonal pinch" mode near 1469 cm"\ All spectra were 
analyzed using Peakfit® software (Seasolve Software, Inc.), assuming a Voigt 
lineshape. Line positions, widths and shapes were allowed to vary to obtain the best 
possible overall fit, and the peak areas obtained were used to find the relative 
intensities of the lines analyzed. 
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